Probing the interaction between Rydberg-dressed atoms through interference 
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We study the dynamics of an atomic Bose-Einstein condensate in an optical lattice in which 
the electronic groundstate of each atom is weakly coupled to a highly excited Rydberg state by 
a far off-resonant laser. This dressing induces a switchable effective soft-core interaction between 
groundstate atoms which, in the lattice, gives rise to on-site as well as long-range interaction terms. 
Upon switching on the dressing laser the Bose-Einstein condensate undergoes a nontrivial collapse 
and revival dynamics which can be observed in the interference pattern that is created after a release 
of the atoms from the optical lattice. This interference signal strongly depends on the strength and 
the duration of the dressing laser pulse and can be used to probe and characterize the effective 
interaction between Rydberg-dressed atoms. 



I. INTRODUCTION 

Ultracold atoms trapped in optical lattices provide a 
highly versatile toolbox for the exploration of the statics 
and dynamics of many-body quantum systems [lj. One 
important experimental demonstration is the paradig- 
matic Bose-Hubbard (BH) model [2J, in which bosonic 
atoms undergo the superfluid-Mott insulator quantum 
phase transition. This quantum phase transition is 
driven by the competition of atomic tunneling between 
lattice sites and short-range (on-site) two-body interac- 
tions. Recently, there is a growing interest in the study 
of extended BH models, in which the two-body interac- 
tion is long-range in the sense that it extends over several 
lattice sites. Extended BH models exhibit a host of quan- 
tum phases, amongst them supersolid and checkerboard 
phases [3j H] . Many studies employ dipolar atoms [5] or 
polar molecules [5] to establish many-body quantum sys- 
tems with long-range interactions. A review about recent 
progress along this direction can be found in Ref. 

Extraordinarily strong and long-range multipolar in- 
teractions are also present between atoms excited to Ry- 
dberg states. Here the interaction strength can exceed 
that between groundstate atoms by more than ten or- 
ders of magnitude. The coherent excitation of Rydberg 
atoms by lasers has recently been demonstrated in a se- 
ries of experiments (for a review see [7J and references 
therein) and many theoretical and experimental groups 
have studied the strongly correlated many-body dynam- 
ics of Rydberg gases [8HT2"] and Rydberg atoms confined 
to optical lattices [I3H2T] . 

The long-range interaction between Rydberg atoms is 
usually not of direct use in the context of extended BH 
models as the atomic motion due to coherent tunneling 
between the sites usually takes place on a time scale that 
is a hundred to a thousand times longer than the life- 
time of the Rydberg states. This limitation can be over- 
come by using Rydberg-dressing, i.e. a weak admixture 
of a Rydberg state to the electonic groundstate, rather 
than directly exciting Rydberg states resonantly. Such 
Rydberg-dressed groundstate atoms exhibit long-range 
interactions that are comparable with typical ultracold 
energyscales [22 24 . A number of recent studies have ad- 
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FIG. 1. (a) Schematics of the two-photon coupling between 
groundstate and Rydberg state. The laser which couples the 
state \g) to the intermediate |P)-state is far detuned from res- 
onance such that the groundstate \g) and the Rydberg \nS) 
state form an effective two-level system. Here fl is the effective 
Rabi frequency and A the (two-photon) detuning of the exci- 
tation laser frequencies with respect to the atomic transition, 
(b) Effective interaction potential between dressed atoms and 
its length scales in relation to the optical lattice. At short dis- 
tances (R <C -Rc) the interaction potential is constant and at 
large distances (fl> R c ) it is of van-der Waals type. (c,d,e) 
Envisioned experimental sequence and timings. Firstly (c), 
a BEC is prepared in an optical lattice. Secondly (d), an 
off-resonant laser coupling the groundstate to the Rydberg 
state is applied for a certain time. After the excitation laser 
pulse, the optical lattice is switched off immediately and the 
dressed groundstate atoms expand freely. Finally (e), inter- 
ference patterns are recorded by taking absorption images of 
the expanded atomic cloud. 



dressed the dynamics in Rydberg-dressed atomic gases, 
such as the formation of supersolids [23 26 , solitons [S7J 
and collective many-body excitations [22) as well as ex- 
citation transport [28j . 

Here we consider the far off-resonant laser dressing 
of the electronic groundstate with a Rydberg nS state, 
where n is the principal quantum number and S cor- 
responds to the electronic angular momentum quantum 
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number I — 0. The laser excitation is achieved by a two- 
photon process, as illustrated in Fig. [TJa,. In the regime 
where the effective Rabi frequency f2 is much smaller than 
the detuning A the laser admixes merely a fraction of the 
Rydberg nS state to the groundstate. As shown in Refs. 
[22H24] this induces an effective interaction potential be- 
tween atoms which is sketched in Fig. [TJd. The hallmark 
of this potential is the characteristic change of its behav- 
ior in the vicinity of the distance i? c , which - depend- 
ing on laser detuning and the Rydberg state - typically 
ranges from 1/im to a few /im: At large inter-atomic dis- 
tance R ^> R c , the potential decreases proportional to 
1/i? 6 , originating from the van der Waals (vdW) inter- 
action between Rydberg atoms in nS state. At short 
distances R <C R C1 the interaction potential levels off 
to a constant value. This soft-core behavior is a conse- 
quence of the Rydberg blockade effect [29] , which inhibits 
the simultaneous excitation of nearby atoms to Rydberg 
states. 

In this work we will investigate the effect of the Ryd- 
berg induced interaction on a lattice gas of groundstate 
atoms that is prepared in a superfluid state. In particu- 
lar we are interested in studying the collapse and revival 
dynamics of the interference pattern that emerges upon 
the release of the atoms from the lattice. It has been 
demonstrated that collapses and revivals of a superfluid 
state released from an optical lattice provide important 
information on two-body interactions (301 131) . For exam- 
ple, complete and periodic revivals are expected if atoms 
interact through short-range interactions (30] . In this 
work, we find that the peculiar shape of the interaction 
potential between Rydberg-dressed atoms gives rise to 
characteristic features in the collapse and revival dynam- 
ics of the interference pattern which originate from the 
competition of short-range and long-range parts of the 
interaction. Interference experiments therefore provide a 
way to probe the presence and to characterize the effec- 
tive interaction between Rydberg-dressed atoms. 

The paper is organized as follows. In Sec. [nj we de- 
rive an effective extended BH model that governs the 
dynamics of Rydberg dressed atoms in an optical lattice. 
In Sec. |III| we study the evolution of a superfluid state 
under the influence of the dressing, using the superfluid 
order parameter. In Sec. IV we analyze the interference 
pattern of dressed atoms released from a one- and two- 
dimensional optical lattice. We conclude and provide ex- 
perimental parameters in Sec. [V] 



II. EFFECTIVE HAMILTONIAN 

The system we consider here consists of N bosonic 
atoms distributed over L D lattice sites, where D is the 
spatial dimension. In our treatment, the Hamiltonian of 
the dressed groundstate atoms is obtained in two steps. 
First, we calculate the Born-Oppenheimer many-body in- 
teraction potential by a perturbative diagonalization of 
the electronic Hamiltonian. With this interaction poten- 



tial, we then derive an extended BH Hamiltonian that 
governs the external dynamics of the dressed groundstate 
atoms in the optical lattice. 

The internal level structure of the atoms is modeled 
by two states, the electronic groundstate \g) and the Ry- 
dberg \nS) state. These two electronic states are cou- 
pled by a two-photon transition with an effective Rabi 
frequency fl as shown in Fig. [IJi. The Hamiltonian for 
the internal (electronic) degrees of freedom of the atomic 
ensemble is given by (using h = 1 and the rotating-wave- 
approximation) 



V v dw( f rj,rk)\nS)j(nS\ ® \nS) k (nS\ 



where H = Q(\g)j(nS\ +h.c.) + A\nS)j(nS\ is the single 
atom Hamiltonian. Here Kdw( r j, r fc) = Ce/\ r j — r-fc| 6 is 
the vdW interaction that is present between two Ryd- 
berg atoms at positions Vj and . Typical values for the 
dispersion coefficient C§ in the case of rubidium-87 can 
be found in [33] . We do not consider the two-body 
contact interaction between groundstate atoms. This in- 
teraction can be switched off for instance by a Feshbach 
resonance [53] . 

In our system, the electronic dynamics (Rydberg ex- 
citation and van-der- Waals interaction) takes place on 
a timcscalc (~ /is) that is orders of magnitude faster 
than that of the external motion of ultracold atoms (~ 
ms). The huge difference in these time scales permits the 
use of the Born-Oppenheimer approximation for treating 
the electronic parts of the Hamiltonian. Rydberg dress- 
ing implies a far off- resonant excitation, i.e., |fi/A| <C 1 
which allows to adiabatically eliminate the atomic Ryd- 
berg states. Using A > we obtain in fourth order per- 
turbation theory (in the small parameter O/A) a Born- 
Oppenheimer surface which can be written as a constant 
plus a sum of two-body interactions [23] HI] > 
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with 



V(r j) r k ) = 



>R? 



Here g = 2ft 4 /A 3 and R c = (C6/2A) 1 / 6 represents the 
characteristic length scale of the soft-core interaction po- 
tential sketched in Fig. In the following, we will ne- 
glect the constant terms in Eq. (JlJ which correspond to 
the second and fourth order light shift. 

Let us now turn to the discussion of the external (mo- 
tional) degrees of freedom. The Rydberg dressed atoms 
are trapped in a Z)-dimensional optical lattice potential 
Vl (r) (lattice spacing d) which is experimentally created 
by an optical standing wave [30 . For D < 3 we assume 
that there is a tight transverse confining potential with 
respect to which the atoms are in the motional ground- 
state. We furthermore neglect differential dephasing of 
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the electronic groundstate with respect to Rydberg state 
which would be caused by different effective trapping fre- 
quencies due to varying AC-polarizabilities of the two 
states. Such dephasing of the motional degrees of free- 
dom can be minimized by choosing a particular frequency 
of the standing wave light field that produces the optical 
potential in conjunction with an appropriate choice of 
the Rydberg state (TI5]. With these approximations the 
second quantized Hamiltonian for the external motion of 
the dressed atoms reads 

r h 2 

H= / dr*t( r )[ V 2 + Vi(r)]*(r) 

J 2m 

drdr / * t (r)* t (r / )F(r,r / )*(r')*(r) (2) 

where \&(r) is the bosonic field operator of the dressed 
atoms. We expand the field operator \l/(r) in terms of 
Wannier states w(r) in the lowest Bloch band of Vl(i") 
such that ty(r) = J2j^j w ( r ~ r j) where bj is the anni- 
hilation operator of a dressed atom at site j and Vj is 
the spatial position around which the Wannier states are 
centered. We consider the case of a deep optical lattice 
where these Wannier states are strongly localized on a 
single lattice site, i.e. the variance of the Wannier wave 
function is much smaller than the lattice spacing d. Here 
the Hamiltonian becomes that of an extended BH model 



U,k) 3 



(3,k) 



with (j, k) indicating the summation over neighboring 
lattice sites and rij — bjbj. The constants J and g refer 
to the tunneling matrix element and the on-site short- 
range interaction, respectively, while 

gR e c 

1[m) |m| 6 d 6 + i?6 

characterizes the long-range interaction among atoms lo- 
cated in distant sites j and k. Note that in general the 
index j is a vector for D > 1. For example, j = {j x ,jy,jz) 
in a three-dimensional optical lattice. 



III. COLLAPSE AND REVIVAL OF THE 
SUPERFLUID ORDER PARAMETER 



In this section, we study the effect that the Rydberg 
dressing has on the coherence properties of the atoms. 
The situation we have in mind is depicted in Fig. [l]>e. 
At first (Fig. [TJ;) one prepares - in the absence of the 
dressing laser - a superfluid state which is the ground- 
state of Hamiltonian (|3| with 5 = 0. In the second step 
(Fig. [TJi) one applies the dressing laser for a time t. Here 
we consider a regime in which the parameters are chosen 
such that |J| <C |g|,|7(m)|, i.e. the atomic tunneling 
can be neglected during the time the dressing laser is 
switched on. The system then evolves approximately un- 
der the Hamiltonian ^ with J = 0. We are interested in 



the coherence properties of the superfluid state after this 
procedure. To this end we will study the superfluid order 
parameter in this section, while in the next section the 
probing of the coherence properties via an interference 
experiment (Fig. [l^) is discussed. 

The initial superfluid state of the atoms is given in 
terms of many-body Fock states, 



|SF) = 



1 



|0,0, 



(4) 



where |0, 0, ■ • • ) is the product of vacuum states of sites 
j = 1,2,.... We work in the limit of a large atom number 
N where we can replace the exact superfluid state by 



isF)«nK-> 



(5) 



with \ctj) being a coherent state of atoms in a single lat- 
tice site defined as, 



E 



(6) 



We will use this representation as initial state for the 
following calculations and will moreover exclusively con- 
sider a homogeneous system where the average atom 
number in each lattice site is constant, i.e., otj = a with 
H 2 = N/L. 

When the Rydberg excitation laser is applied, 
the superfluid state evolves according to |SF(£)) = 
exp[— iHt]\SF(t — 0)). We will now discuss the effect 
of the dressing on the superfluid order parameter defined 
as 



(-D) 



(*) = (SF(t)|6 3 -|SF(t)) = (bj) = aexp [\a\ 2 F D (t)] . 



We consider D-dimensional optical lattices with linear 
(D = 1), square (D = 2) and cubic (D — 3) geometry. 
We will focus our analysis on the "phase factor" Fu(t) 
as it encodes the full dynamics of the order parameter 
4>( D \t) and has the advantage that it is independent of 
the particle density. For the following discussion we will 
furthermore consider a parameter regime in which the 
interactions between distant atoms, i.e. atoms whose po- 
sition labels obey \j — k\ > 2, is negligible. This is the 
situation depicted in Fig. [T]d. However, our results can 
straightforwardly be generalized. Let us now study Fr>(t) 
for lattices with different spatial dimension D: 
(i) One- dimensional lattice. Taking into account the 
on-site, nearest-neighbor 7(1) and next-nearest-neighbor 
7(2) interaction we find 

Fi(t) = -[5 - c- l9t - 2c~ l ^ t - 2e-^ 2 )*]. (7) 

In general we expect the two-body interactions to result 
in a collapse and revival of the superfluid order [SUl [3S] . 
The collapse dynamics can be obtained by analyzing the 
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real part of F\(i) in the vicinity of t = 0. A Taylor 
expansion up to second order in t yields 



Re[Fi(t)]«- 



+ 7 2 (l) + 7 2 (2) 



t 2 = - 



(8) 



where t± is the collapse time, given as 



1 



n 



v / ff 2 + 27 2 (1)+ 27 2 (2) ' 



This expression clearly indicates that the collapse is 
faster the stronger the two-body interaction. The pos- 
sibility of revivals of the initial state is strongly deter- 
mined by the long-range part of the interaction. Com- 
plete revivals occur when all phase factors appearing in 
exponential terms of F\ (t) are simultaneously multiples 
of 2tt. As the ratios between the various possible interac- 
tion strengths are in general irrational, complete revivals 
are unlikely and one can rather expect partial revivals at 
a finite time. Before discussing this aspect in more de- 
tail let us briefly provide the results for two- and three- 
dimensional lattices. 

(ii) Two-dimensional lattice. In two dimensions the trun- 
cation of the interaction to atoms whose indices obey 
\j — k\ < 2 means that we have to take into account 
at most next-next-nearest neighbor interactions. This 
yields 



F 2 (t) = -[13 



-igt 



4e 



-i 7 (l)t 



4e 



4e 



and a collapse time 

72 



^. 9 2 + 4 7 2 (1) + 4 7 2 (V2) + 4 7 2 (2) 

(iii) Three-dimensional lattice. In this case one obtains 

F 3 (t) = -[33 - e~ l3t - &e- 1 ^ 1 - 12e~ i7(v/2)t 
_ 8e -n(V3)t_ 6e -n(2)t] ! 

with a corresponding collapse time 

1 

yV 2 + 6 7 2 (1) + 12 7 2 (V2) + 8 7 2 (V3) + 6 7 2 (2) 

Generally the collapse time becomes shorter for higher 
dimensional systems, since the number of contributing 
interaction terms grows with increasing dimensions (as 
a consequence of the increasing coordination number). 
For the purpose of illustration let us now analyze the 
case D = 2 in more detail. In particular, we intend 
to study F2(t) as a function of the ratio between the 
lattice spacing d and the length scale R c which demar- 
cates the soft core of the effective interaction potential 
among Rydberg-dressed atoms. In Fig. [2] we present 
data for d/R c = {0.7,1.0,1.5}. The smaller d/R c the 
more terms of the long-range tail of the interaction actu- 
ally contribute to the denominator of t<x- This leads to 
a more rapid collapse of the order parameter the smaller 
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FIG. 2. (a) Dynamical evolution of F2(t) for various values 
of d/R c . (b, c) Dotted curves show calculations merely tak- 
ing into account the on-site and nearest-neighbor interaction. 
Solid lines show calculations including all interactions up to 
next-next-nearest neighbors. 



d/R c becomes, which is clearly shown by the data. In 
fact, all dimensions have in common that for d/R c > 1, 
the behavior of the function Fjj (t) is determined mainly 
by the on-site and nearest-neighbor interaction. Com- 
pared to those the other interaction terms are negligible 
which is owed to the fast drop of the soft-core interaction 
potential at large inter-site separations. To demonstrate 
this explicitly, we calculate F%(t) by neglecting interac- 
tion beyond nearest neighbors, i.e., 7 (m) = for m > 1. 
The results and a comparison with the full calculation 
are shown in Fig. and c. For d/R c = 0.7 the two cal- 
culations deviate already beyond gt = 2. However, for 
d/R c = 1.5 agreement for much longer times (here until 
gt = 30) is achieved. Thus, for sufficiently short times 
the contribution of the long-range interaction terms to 
the dynamical evolution of i 7 ^) is strongly suppressed. 



IV. INTERFERENCE PATTERN 

The dynamical evolution of the superfluid order can 
be made visible from an interference experiment as was 
shown in Ref. 30 . To create the interference pattern 
one simultaneously switches off the optical lattice and the 
Rydberg excitation laser. After a time-of-flight period, 
the interference pattern in the spatial density distribu- 
tion (see Fig. [!£) is formed and made visible by taking 
an absorption image of the atomic cloud. This inter- 
ference pattern is characterized by the quasi-momentum 
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distribution [13 137] 

L 

S(k) = J2^ irm - rn) (blb n ). (9) 

ran 

Here k is the quasi-momentum with components = 
V£ x 2ir/ (dL) where = 1, 2, • • • , L and, e.g., £ = x,y, z 
when D = 3. For convenience we furthermore define the 
dimensionless quantity k = c?k. 

Let us now study the effect of the Rydberg-dressing on 
5(k). We will consider here merely on-site and nearest- 
neighbor interactions of the extended BH model. This 
approximation is well justified provided that d/R c > 1, 
as shown in the previous section. We furthermore will 
only focus on one- and two-dimensional lattices. In or- 
der to calculate S(k), we first evaluate the matrix el- 
ements of the single particle density matrix, (b^bk) = 

<SF(i)|&t& fc |SF(t)) which are given in Appendix A Using 
this data we obtain the following results: 

(i) One- dimensional lattice. In order to be able to distin- 
guish contributions of the on-site and nearest-neighbor 
interaction, we start with a situation where only the on- 
site interaction is present. This corresponds to a sce- 
nario that was experimentally studied in Ref. |30j and in 
case of the Rydberg-dressed atoms this is realized when 
d/R c 3> 1. The quasi-momentum distribution S(k x ) 
changes as a function of time as shown in Fig. [3^. Two 
main features can be identified. First, at t = 0, S(k x ) is 
centered around k x = 0. As the interaction time grows 
S(k x ) collapses to flat distribution and at an even later 
time (gt = 2tt) it revives to its original value. The sec- 
ond feature is that the periodicity of S(k x ) does not de- 
pend on k x (although the amplitude does). This is a 
direct consequence of the fact that the correlation func- 
tion is factorizable, i.e., = for all j and k 
[55] , The time-dependence thus becomes a global factor 
of S(k x ). 

This changes in the presence of the nearest-neighbor 
interaction. Due to a nonzero 7(1), (bjbk) is no longer 
factorizable but depends on the relative distance be- 
tween respective sites (Appendix |A| . This gives rise to 
strikingly different features in the quasi-momentum dis- 
tribution. An example of this is shown in Fig. ^jp for 
d/Rc = 1.5 where 7(1) « 0.08g. The noteworthy feature 
here is that the time-dependence of S(k x ) depends on 
the quasi-momentum k x . For example, at gt — 2n, we 
find that S(k x ) revives partially along k x = (Fig. [3j;), 
while along k x = ir it shows a clearly different oscillation 
pattern (Fig.|3ji). 

(ii) Two-dimensional lattice. Let us now turn to the 
discussion of the interference pattern that emerges in 
the case in which the atoms are released from a two- 
dimensional optical lattice. In order to have a reference 
we again calculate at first the interference patterns for 
vanishing nearest-neighbor interaction. The main fea- 
tures in the interference pattern can be captured by ana- 
lyzing time slices taken at gt = 0, n/2 and 2ir. This data 
is shown in Fig. [4^,-c. Like in the one-dimensional case 
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FIG. 3. Quasi-momentum distribution of N = 82 atoms re- 
leased from a one-dimensional optical lattice with L — 41 
sites, (a) On-site interaction only, i.e. y(m) = (correspond- 
ing to d/R c ^> 1), (b) On-site and nearest-neighbor interac- 
tion for d/R c = 1.5 (7(1) ~ O.O85). (c,d) Cuts through S(k x ) 
at different values of the quasi-momentum. The panels show 
cuts along k x = (c) and k x = it (d). The solid/dashed 
curves correspond to the parameters of panel (a) /(b). 
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FIG. 4. Natural logarithm of the quasi-momentum distribu- 
tion S(k) of atoms released from a two-dimensional 15 x 15- 
lattice and an average particle number in each sites of |a| 2 = 
2. The figures show snapshots taken different times: gt = 
[(a) and (d)], gt — n/2 [(b) and (e)] and gt = 2n [(c) and 
(f)]. Panels (a)-(c) show data that has been obtained by tak- 
ing only on-site interactions into account. The data shown in 
panels (d)-(f) includes also nearest-neighbor interactions. 



the time-evolution is periodic showing a first revival at 
gt = 2-k and the expected periodicity [3D]. In Fig. [4ji-f, 
we present the interference pattern that is obtained in 
case of a finite nearest-neighbor interaction (d/R c = 1.45 
and 7(1) ~ O.lg) using the same time slices. We find that 
the height of the peaks of 5(k) in Fig. [4^ is lowered com- 
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pared to Fig. [4j;. The most striking feature is, however, 
the appearance of additional interference maxima in be- 
tween the ones that occurred previously in the case of 
pure on-site interactions. These peaks occur because the 
nearest-neighbor interaction actually affects (indirectly) 
the phase relation between atoms that are separated by 
a distance y/2 d (next-nearest neighbors) as shown in Ap- 
pendix [XJ 

V. CONCLUSIONS AND OUTLOOK 

In conclusion, we have shown that Rydberg-dressed 
atoms in an optical lattice give rise to an extended BH 
model. In particular, we have studied the dynamical evo- 
lution of an atomic superfluid under the influence of Ry- 
dberg dressing. The emerging long-range interactions re- 
sult in a rapid collapse of the superfluid order parameter 
and in general allow only for partial revivals. Moreover, 
we have demonstrated that interference experiments can 
directly reveal the interaction between Rydberg-dressed 
atoms. The hallmark here is the emergence of additional 
interference maxima. In addition, the time-dependence 
of the quasi-momentum distribution is no longer a global 
factor. This leads to a time-dependence which strongly 
depends on value of the quasi-momentum. 

We expect that our predictions can be probed in cur- 
rent experiments using Rydberg atoms in lattices, e.g. 
the one presented in Ref. [T3] by Viteau et al. In this 
particular experiment Rydberg D-states of rubidium-87 
with principal quantum numbers in the range n = 55 — 80 
are excited and the lattice spacing can be tuned within 
0.42 [im < d < 13 /*m. If one instead excites a Rydberg 
S-state with n = 60 and uses a laser with the parameters 
A = 2-7T x 15 MHz and £1 = 2ir x 0.5 MHz one obtains 
R c fa 4 /im as the characteristic length scale of the soft- 
core potential. This leaves sufficient freedom for tuning 
the ratio d/R c . The typical revival time of the inter- 
ference signal evaluates to to ps 27 ms. This has to be 
compared with the effective lifetime of dressed ground- 



state atoms, r G ff = |A/fi| 2 x TR y d where TR y( j is the bare 
lifetime of the Rydberg state. Using the above parame- 
ters, we obtain T e g = 225 ms, which leads to a coherence 
time that in principle should permit the observation of 
the predicted dynamics. 
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Appendix A: Correlation functions 

Using the initial state Eq. ^ and taking into ac- 
count of on-site and nearest-neighbor interactions, (bjb k ) 
is calculated analytically. The matrix elements required 
for calculating the quasi-momentum distribution are pre- 
sented below. 

If j = k, we get (bjbj) = \a\ 2 . 

If \j -k\ = l, 

If \j - k\ = V2, 

(b)b k ) = \ a \ 2 e~ 2]a]2[1 ~ COS9t] ~ i(D ~ 1)]a f [1 ~ cos ^ {1)t] . 
If \j -k\=2, 

(b)b k ) = | Q ,| 2 e - 2 l"| 2 [ 1 - cos 9 t ]-2(2£'-l)|Q| 2 [l-cos7(l)t]^ 

If \j -k\>2, 

(btb k ) = | a | 2 e" 2|Q|2[1 ~ COS9tl e" 4I)|a|2[1 " COS7(:L) * ] . 
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